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Abstract 

We present a measurement of the branching ratio of the CP violating decay -RTl— > 
7r"'"7r~ performed by the KLOE experiment at the (p factory DA<I>NE. We use 328 
pb~^ of data collected in 2001 and 2002, corresponding to ~ 150 million tagged 
Kl mesons. We find BR{Kl Tr+vr") = (1.963 ± 0.012stat. ± 0.017syst.) x 10"^. 
This branching ratio measurement is fully inclusive of final-state radiation. Using 

the above result, we determine the modulus of the amplitude ratio \rj^ | to be 

(2.219 ± 0.013) X 10-3 and |e| to be (2.216 ± 0.013) x lO'^. 



CP violation was discovered in 1964 through the observation of the decay 
TT+TT" [1]. The value of JiYi^Ki tt+tt-) is known today with high accuracy from 
the results of many experiments, but a recent measurement by KTeV [2] is in dis- 
agreement with the value reported by the PDG [3]. In the Standard Model, CP 
violation is naturally accommodated by a phase in the quark mixing matrix [4,5]. 
BR(Jri— ^TT+TT-), together with the well known values of BR(ii'5— S'Tt+tt-), tksi 
and tkl, determines the modulus of the amplitude ratio |?7+-| = 
^Jt{Kl —> 7i~^7i~) /T{Ks vr+TT-), which is related to the CP violation param- 
eters e and e' by = e + e' ~ e. The value of |e| can thus be obtained and 
compared with the Standard Model calculation. In this letter, we present a mea- 
surement of the branching ratio of the decay Kl tt^tt^ performed by the KLOE 
experiment at the factory DA$NE. 

The KLOE detector consists of a large, cylindrical drift chamber (DC), sur- 
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rounded by a lead/scintillating- fiber electromagnetic calorimeter (EMC). A su- 
perconducting coil around the calorimeter provides a 0.52 T field. The drift 
chamber [6] is 4 m in diameter and 3.3 m long. The momentum resolution is 
c"px/P-L ~ 0.4%. Two-track vertices are reconstructed with a spatial resolution 
of ~ 3 mm. The calorimeter [7] is divided into a barrel and two endcaps. It 
covers 98% of the solid angle. Cells close in time and space are grouped into 
calorimeter clusters. The energy and time resolutions for photons of energy E 
are as/E = 5.7%/ ^Je (GeV) and at = 57 ps/^E (GeV) © 100 ps, respectively. 
For this analysis, events triggered [8] only by calorimeter signals are used. Two 
energy deposits above threshold {E > 50 MeV for the barrel and > 150 MeV 
for the endcaps) are required. 

For the present measurement, we use a subset of the data collected by KLOE 
during the years 2001 and 2002 that satisfies basic quality criteria [9]. The data 
used corresponds to an integrated luminosity of ~328 pb~^. Each run used in 
the analysis is simulated with the KLOE Monte Carlo program, GEANFI [10], 
using values of relevant machine parameters such as -^s and as determined 
from data. Because of the small beam-crossing angle, mesons are produced 
at DA$NE with a momentum of about of 12 MeV/c, toward the ring's center. 
Machine background obtained from data is superimposed on Monte Carlo events 
on a run-by-run basis. The number of events simulated for each run in the data set 
is equivalent to that expected on the basis of the run luminosity. For this analysis, 
we use a Monte Carlo sample consisting of (p KsK^ events in which the Kg 
and Kl decay in accordance with their natural branching ratios. The effects of 
initial- and final-state radiation are included in the simulation. The treatment 
of final-state radiation in Ks and decays follows the method discussed in 
Ref. 11. In particular, the >7r+7r~(7) event generator includes the amplitudes 
for inner bremsstrahlung and direct emission. Good runs, and the corresponding 
Monte Carlo events, are organized into 14 periods. The branching ratio analysis 
is performed independently for each of these periods. The average over the 14 
periods gives the final result. 

At a factory, neutral kaons are produced through (p — > KsKl decays. A pure 
sample of K^s can be selected by identification of Ks — tt+tt" decays (tagging). 
The tagging efficiency depends slightly on the fate of the Kl. it is different for 
events in which the Kl decays to each channel, interacts in the calorimeter, or 
escapes the detector. KLOE has already used the tag method to measure the 
dominant Kl branching ratios and the Kl lifetime [12]. The tagging criteria de- 
scribed in Ref. 12 where chosen to minimize the difference in the tagging efficiency 
among the various decay modes. In the present analysis, in order to increase the 
statistics, we employ a tagging algorithm with an efficiency of about a factor of 
seven larger than that used in Ref. 12. 

A precise measurement of BR(A'l tt+tt^) is obtained here by measuring the 
ratio R = BR{Kl 7r+7r^)/BR(A'i vr^/i^i^) since the values of of the tagging 
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efficiency for Kl^tt^tt and Kl-^ti^v are similar. The value of BR(i^/^ — > tt+tt ) 
is finally obtained using the value of BR(ii'2. — * n'^fi^h') from Ref. 12. 

For J^^. tagging, we select Ks tt+tt" decays by requiring two tracks of opposite 
curvature from a vertex in a cylindrical FV with r^y < 10 cm and \z\ < 20 cm, 
centered on the collision region as determined for each run using Bhabha events. 
We also require that the two tracks give |m(7r"'"7r~) — mKo\ < 5 MeV/c^ and that 
||p_,_ + p_| — Pks\ < 10 MeV/c, with p^s calculated from the kinematics of the 
(j) — > KsKl decay. The momentum, Pkli is obtained from the Ks direction 
and The Kl line of flight {tagging line) is then constructed from the Kl 
momentum, y>Kl = P-ftTs' ^"^^ position of the production vertex, x,^. 

The overall value of the tagging efficiency is about 66%, and that of the ratio of 
the tagging efficiencies for events in which Kl^ti^tt~ and Kl^tt^v, averaged 
over all periods, is 1.012 ± 0.001 as obtained from Monte Carlo. 

To reconstruct the decay vertex, all relevant tracks in the chamber, after 
removal of those from the Ks decay and their descendants, are extrapolated to 
their points of closest approach to the tagging line. For each track candidate, we 
evaluate the point of closest approach to the tagging hne, Xc, and the distance of 
closest approach, dc- The momentum pc of the track at Xc and the extrapolation 
length, /c, are also computed. Tracks satisfying dc < ar^y + fc, with a = 0.03 and 
6 = 3 cm, and —20 cm < 1^ < 25 cm are accepted as K^ decay products, r^y 
being the distance of Xc from the origin in the transverse plane. For each sign 
of charge we consider as associated to the Kl decay the track with the smallest 
value of dc- We attempt to find a vertex using these two tracks. Candidate Kl 
vertices must be within the fiducial volume defined by 30 < r^^y < 150 cm and 
\z\ < 120 cm. 

The tracking efficiency has been evaluated by Monte Carlo with corrections ob- 
tained from data control samples [9]. The corrections range between 0.99 and 
1.03 depending on the i^L-decay channel and on the run period. We determine 
the conditional tracking efficiency as the ratio of the number of events in which 
there are two i^x-decay tracks of opposite sign and the number of events in 
which there is a single -R'^.-decay track. We use the same method for both data 
and Monte Carlo, and correct the efficiency obtained from the simulation with 
the data-Monte Carlo ratio of the conditional tracking efficiencies from control 
samples. The corrections are evaluated as functions of the track momentum using 
Kl-^h^7i^it^ and Kl—^ticv control samples. To select Kl^tx^ti^tx^ events, we 
require at least one K^ decay track and two photons from decay. To reject 
background, we apply cuts on the two-photon invariant mass and on the photon 
times of flight. The momenta of decay tracks are calculated with a resolution of 
about 10 MeV/c using the two photons from the 7r° decay and the momentum of 
the other track. In order to obtain the corrections for higher momentum tracks 
(above ~150 MeV/c), Ki^neu decays are used. To select Kl^txcu decays, we 
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identify electrons by time of flight with a purity of about 95%. The background 
is mostly due to decays. The momentum of the second track is evaluated 
from the missing energy with a resolution of about 30 MeV/c. The value of the 
tracking efficiency is about 61% for 7r"'"7r~ decays and 51% for ir^fi^u decays. 
Once two tracks have been found, the vertex efficiency is about 97%. We find 
values for data and Monte Carlo to be in agreement within 0.1%. 

The number of —>■ n^fi^i'l'j) events is obtained by fitting the distribution of 
^TTA* = IPmissI ~ -S'miss with a linear combination of Monte Carlo distributions as 
in Ref. 12. Att^^ is obtained from the smaller absolute value of the two possible 
values of |pmiss| — -E'miss, where Pmiss and -Emiss are the missing momentum and 
missing energy at the Kl vertex, evaluated by assuming the two tracks to be a 
TTyU or /ivr pair. 

The best variable for selecting Kl — > 7r+7r~(7) decays is E'^^^^ + |pmissP where 
i?miss in this case is the missing energy in the hypothesis of the tt+tt^ 
decay. The result of the fit to the y^E'^j^g + |pmissP distribution with a linear 
combination of Monte Carlo distributions is shown in Fig. 1, left for a single run 
period. The signal region is shown expanded in Fig. 1, right. The sum over the 

^5000f 1 ^ F 1 
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Fig. 1. Distribution of ^J E'^^^^^ + |pmissP for a single run period, with fit to Monte Carlo 
distributions for different decay channels: full fit range (left panel); signal region (right 
panel). 

14 periods of the numbers of signal and normalization events obtained from the 
fits, and the average values over the 14 periods for the efficiency ratios are given 
in Table 1. 

The systematic errors due to the limited knowledge of the corrections, data-Monte 
Carlo discrepancies, and instabilities in the signal counting have been studied as 
described in the following. 
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number of ir~^Tr events 


45267 ± 255 


number of tt^u events 


5243723 ± 2805 


^tagging ('^"'"'^ )/£tagging ('^M'^) 


1.0119 ± 0.0008 


etracking(7r"'"7r~)/ etracking (TT/il^) 


1.1726 ± 0.0011 



Table 1 

Analysis summary: number of signal and normalization events obtained from the fits, 
and the average values of the efficiency ratios. 

The systematic error on the determination of the ratio of the tagging efficiencies 
is evaluated from the the stability of the results obtained with different tagging 
criteria. Since the variations in the tagging efficiency are mostly due to the de- 
pendence of the calorimeter trigger on the Kl decay channel, the analysis has 
been repeated with the additional requirement that the event trigger be due to 
the Ks decay products alone {self-triggering event). Another contribution to the 
systematic uncertainty on the tagging efficiency is due to the dependence of the 
reconstruction efficiency of the pions from Ks ir^ir" on the presence of other 
tracks in the drift chamber. This effect is reduced for events in which the pion 
directions are almost orthogonal to the Ks line of flight. The analysis has been 
repeated with the requirement that the difference between the Ks pion momenta 
be less than 60 MeV. The largest fractional difference in the results obtained 
using different tagging criteria is 32 x 10~^. This value is taken as the systematic 
uncertainty on the tagging efficiency. 

The uncertainty on the tracking efficiency is dominated by the control sample 
statistics and by the variation of the results observed by using different criteria 
to identify tracks from Kl decays. The systematic error due to the possible bias 
introduced in the selection of the control sample has been studied by varying 
the values of the cuts made on dc and Ic when associating tracks to K^ vertices. 
These changes result in a variation of the tracking efficiency by about ±20%. 
The corresponding fractional change in i? is 37 x 10~^. This value is taken as the 
systematic uncertainty on the tracking efficiency. 

The systematic error due to the reliability of the Monte Carlo distributions in the 
variables |pmiss| — -E'miss and ^JE^~+[p^^^^ has been evaluated by studying the 

sensitivity of R to the momentum resolution. The fits to the ^ -^miss + I Pmiss I ^ and 
the I Pmiss I — -E'miss distributions have been repeated using Monte Carlo samples 
in which the momentum resolution has been varied according to the uncertainty 
determined using high purity control samples of the dominant Ki decays. A 
fractional systematic uncertainty of 55 x 10~^ is obtained. 

The flnal value for BK^Kl — > tt'^tt") has been calculated as the average over 
the 14 periods of data taking. Good stability of the results is observed after 
application of the period-dependent tracking efficiency corrections. The value 
of the fit corresponds to a confldence level of 77%; see Fig. 2. Including systematic 
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Fig. 2. Values of R for 14 run-periods. The result of a fit to a constant value is also 
shown. 

uncertainties, we obtain: 



BR{Kl TT + TT") 



(0.7275 ± 0.0042stat ± 0.0054syst) x 10" 



where the statistical error includes the contribution from Monte Carlo statistics 
and the statistical uncertainty on the tracking efficiency correction. Using the 
KLOE result BR{Kl 7r/iz/)=0.2698±0.0015 [12] we find: 



BR{Kl ^ vr+vr-) = (1.963 ± 0.012stat. ± 0.017syst.) x 10"^ 

This result is fully inclusive with respect to final-state radiation including both 
the inner bremsstrahlung and the (CP-conserving) direct emission components. 
The result is in good agreement with the measurement from KTeV [2], ^ (1.975 ± 
0.012) X 10~^, and in strong disagreement with the value reported by the PDG, 
(2.090 ± 0.025) X 10-3 [3]. 

The error on the measurement of BR{Kl — »■ vr/iz/) of Ref. 12 is correlated with 
the errors on BR{Kl vrez/), BR{Kl 7r°7r°7r°), BR{Kl 7r+7r-7r°), and the 

^ The contribution of direct emission is claimed to be negligible in the KTeV result. 
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Kl lifetime. This is because, in the treatment of Ref. 12, the branching ratios 
and the hfetime are determined from hfetime-dependent absolute branching ratio 
measurements using the constraint X^BR = 1 — 0.0036, where 0.0036 represents 
the small contribution due to the sum of the branching ratios for Kl — »■ tt+tt", 
Kl —>■ 7r°7r°, and — > 77 listed in the PDG compilation [3]. As a consequence, 
the error on the value for BIi{Ki tt^tt^) reported here is correlated with 
the errors on the branching ratio and lifetime values reported in Ref. 12. 
Additional sources of correlation arise from common systematic uncertainties in 
the branching ratio measurements. The complete correlation matrix is: 

^1 -0.25 -0.56 -0.07 -0.12 0.25 ^ 
1 -0.43 -0.20 0.51 0.33 
1 -0.39 -0.24 -0.21 
1 -0.09 -0.39 
1 0.11 

V 1 ) 

where the columns and the rows refer, in order, to the BR's for decays to TT^e'^i', 
ir'^fi^i', Tr^TT^vr", tt+tt^tt^ and 7r"^7r~, and Tkj^- Because BR(A'l n^n^) is small, 
the KLOE values for the dominant BR's are essentially unchanged when the new 
KLOE value for BK{Kl n'^n') is used in constraining the sum of the Ki BR's 
to unity. 

The measurement of BR(/('i 7i~^7i~) can be used to determine |?7+-| and |e|, 
correcting for the small contribution of e'. Using the measurements of BIi{Ks 
vr"*"-/?") and tkj^ from KLOE [12,13] and the value of tk^ from the PDG [3], and 
subtracting the contribution of the CP-conserving direct-photon emission pro- 
cess Kl —* 7r+7r~7 [14] from the inclusive measurement of BK{Kl vr+vr"), 
we obtain |?7+_| = (2.219 ± 0.013) x 10""^. Finally, using the world average 
measurement of Re(e'/e) = (1.67 ± 0.26) x 10^^ and assuming arge' = arge, 
we obtain |e| = (2.216 ± 0.013) x 10"'^. This result is in disagreement with 
the value |e| = (2.284 ± 0.014) x lO'^ reported in the PDG compilation [3]. 
The value of |e| can be be predicted from the measurement of other observ- 
ables and compared with measurements. For this purpose, we use the prediction 
|e| = (2.875±0.455) x lO'^ obtained by the UTfit collaboration [15,16], where to 
test the mechanism of the CP violation in the Standard Model, the value of |e| has 
been computed from measurements of the CP-conserving observables AM(Pd), 
^M{Bs), Vub, and Vcb- The probability density function of the prediction of |e| is 
shown in Fig. 3. No significant deviation from the Standard Model prediction is 
observed. Notice that, due to the large uncertainties on the computation of the 
hadronic matrix element corresponding to the K-K mixing, the theoretical error 
is much larger than the experimental error. 
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x10"^ 




Fig. 3. Standard Model prediction of |e| from measurements of |Kjb|/|V"cfe|, and AM{Bd) 
and from the limit on AM{Bs) from Ref. 15. 

In summary, a new, precise measurement of the — > vr+vr" branching ratio has 
been performed. The result: 

BR{Kl ^ TT+TT-) = (1.963 ± 0.012stat. ± 0.017syst.) X 10-3 

is in good agreement with that recently obtained by KTeV [2], and in disagree- 
ment with the value reported by the PDG [3]. The impact of this measurement 
on the determination of the apex of the unitarity triangle has been studied and 
compared with determinations from measurements of other observables. 
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